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We demonstrate that spin-charge separation can occur in two dimensions and note its confluence 
with superconductivity, topology, gauge theory, and fault-tolerant quantum computation. We con- 
struct a microscopic Ising-like model and, at a special coupling constant value, find its exact ground 
state as well as neutral spin 1/2 (spinon), spinless charge e (holon), and Z2 vortex (vison) states 
and energies. The fractionalized excitations reflect the topological order of the ground state which 
is evinced by its fourfold degeneracy on the torus - a degeneracy which is unrelated to translational 
or rotational symmetry - and is described by a Z2 gauge theory. A magnetic moment coexists with 
the topological order. Our model is a member of a family of topologically-ordered models, one 
of which is integrable and realizes the toric quantum error correction code but does not conserve 
any component of the spin. We relate our model to a dimer model which could be a spin SU(2) 
symmetric realization of topological order and its concomitant quantum number fractionalization. 



PACS numbers; 75.10.Jm 



I. INTRODUCTION. 



The advent oLnolyacetylene0 and the fractional quan- 
tum Hall efFectErB showed that quantum number frac- 
tionalization is a robust possibility in condensed matter 
physics. The quantum numbers of the low-energy ex- 
citations of these systems are fractions of those of the 
microscopic degrees of freedom, the electrons. There 
are charge e, spin and charge 0, spin i spin-charge 
separated excitations in polyacetylene and other 11? sys- 
tems. The fractional quantum Hall state at filling frac- 
tion v = 1/m has charge e/m, statistics 7r/m excita- 
tions; more exotic possibilities lurk at other v. 
spite a flurry of interest generated by the suggestior 
that spin-charge separation is the mechanism for high- 
temperature superconductivity in the cuprates, it is, at 
present, unclear whether spin-charge separation can oc- 
cur in a 2D jmagnet. There is a set of long-wavelength 
field-theoriesrilil which describe the properties of puta- 
tive fractionalized magnets, but their existence has been 
controversial for want of a concrete microscopic model 
of spin i moments coupled by short-ranged interactions 
in which fractionalization occurs. In this paper, we con- 
struct such a microscopic model of a 2D magnet. We find 
the exact ground state and neutral, spin i (spinon) and 
charge e, spiiiless, (holon) excited eigenstates as well as 
a Z2 vortexllJtael. ThpJxactionalized excitations reflect 
the topological ordertJOo of the ground state w-is 



evinced by its fourfold degeneracy on the toruai3i 
- a degeneracy which is unrelated to translational or 
rotational .symmetry - and is described by a Z2 gauge 
theoryll3c3'Ll. Our construction implies that fractional- 
ization is a reasonable possibility for magnets with short- 
ranged interactions. Our model is a member of a family 
of models, another of which is integrable and realizes the 
toric quantum error correction codetil. The models are 
related to the quantum dimer modelE3 and lie at the con- 
fluence between superconductivity, topology, gauge the- 



ory, and fault-tolerant quantum computation. 

Our purpose here is to show that such microscopic 
models do exist, at least in principle, so we construct 
a model with the aim that it be deep within a phase 
supporting fractionalized excitations, not that it be a 
realistic description of any ppjticular physical system. 
(As we discuss below, KitaevEil, in beautiful work, has 
constructed an exactly soluble model with many of the 
desired properties, but it does not have any conserved 
quantum numbers, so it is not 'fractionalized' in the 
sense of admitting fractional quantum numbers.) How- 
ever, we insist that our model be expressed in terms of 
spin-1/2 electrons, so that it is truly microscopic. Conse- 
quently, our analysis differs in a number ci£-jpii-|«;spects 
from earlier ones which dealt with modelst2rliiliJ'tj which 
are not, strictly speaking, microscopic-electronic models 
or else relied on various assumptionsEja in to order re- 
duce the microscopic models to effective models which 
exhibit fractionalization. We avoid the need for such as- 
sumptions or modifications (however benign they may 
seem) by endowing our model with the following prop- 
erties which distinguish it from other models which have 
been considered in this context: (1) Ising symmetry, (2) 
translational symmetry which is broken by han(i_and (3) 
adiabatic continuability to an integrable modelEll. 

A real magnet will have many additional complica- 
tions, but these are unimportant so long as it shares 
the key feature of our nuid£]__namely topological order. 
In pioneering work, WencjOllj observed that phases of 
matter in two dimensions with fractionalized excitations 
are not characterized by a local order parameter, in con- 
trast to more familiar phases such as crystals. Rather, 
their universal properties are encapsulated by topologi- 
cal quantum numbers, such as their ground state degen- 
eracy on a torus or higher genus surface, over and above 
any degeneracy which is due to broken symmetry. De- 
generacy which is due to topological order persists in the 
presence of local perturbations such as impurities, which 
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break translational and rotational symmetry. (This ob- 
servation will prove important since it guides us to con- 
struct our model on lattices which penalize states which 
would break translational and rotational symmetry on a 
square lattice.) This is completely different from the two- 
fold degeneracy associated with an Ising antiferromagnet, 
which is removed by the application of a small symmetry- 
breaking field at even one point. Topological order is 
well-established theoretically in the fractional quantum 
Hall effecttEl, where it is manifested by the existence of 
excitations with non-trivial braiding statisticscJ. 

Along with spinons and holons- |a_ spin-charge sepa- 
rated state must have Z», vorticesliJo, which have been 
recently dubbed 'visons'a. Topological order implies the 
existence of a gap in the vison spectrum. It is not nec- 
essary for all other excitations to be gapped (see, for 
instance, the construction of Ref. This is analogous 
to the situation in a conventional ordered state such as 
a superconductor, which can have gapless quasiparticles 
if, for instance, it has d-wave symmetry or impurities. 
They do not preclude a stable SC state so long as there 
is a gap to the creation of vortices. Similarly, in our 
topologically-ordered states the existence of a vison gap, 
A„, is necessary to guarantee the existence of distinct 
topological sectors of the Hilbert space on the torus, as 
we will describe below. We compute the vison gap and 
present evidence that the rest of the spectrum is gapped 
(though, we reiterate, this is not a major issue). The 
intcgrable model in the family is fully gapped. 

The concept of topological order is very attractive the- 
oretically because it is precise, but it is sobering to note 
that it has not been possible, to date, to directly mea- 
sure most of the topological quantum numbers - such 
as the braiding statistics - of a fractional quantum Hall 
state. On the other hand, this very feature has generated 
considerable interest in the use of topologically-ordered 
states for quantum computation. The inaccessibility of 
topological degrees of freedom to local probes insulates 
them against many forms of decoherence, the bete noir 
of the quantum computation progrpii. 

This point was made by KitaevEil in a beautiful pa- 
per in which he constructed a concrete model exhibiting 
the requisite topological order and a fault-tolerant quan- 
tum error correctin g co de which could be implemented 
in it (see also Ref. 24). The intcgrable model in our 
family is equivalent to Kitaev's. For our purposes, the 
model of greater physical interest is the one which con- 
serves and exhibits quantum number fractionaliza- 
tion, which is of intrinsic interest and iftight be relevant 
to high-temperature superconductivityQ. It could also 
prove useful for quantum computing since their spin and 
charge quantum numbers allows for the manipulation of 
spinons and holons. Harnessing the otherwise elusive vi- 
sons also becomes a real possibility if the proposed exper- 
iment of Ref. ^ can be implemented. Finally and per- 
haps most importantly, the energy scale associated with 
topological order in a magnet is likely to be an exchange 
constant J ~ lOOOfC. Thus, a magnet with fractional- 



ized excitations has many attractive features as a milieu 
for quantum computation (for another, see Ref. |2^ and 
refs. therein). 



II. THE MODEL. 

Our model has spin ^ degrees of freedom, S^, living on 
the links of a lattice which we specify below. They are 
not gauge fields, but gauge- invariant, physical degrees of 
freedom which happen to be located on the links of the 
lattice (we return to this point later). The Hamiltonian 
is: 

H^ = JiY.g{St) - J2Y.F,Pp + J,Y,Pp (1) 

i P P 

where Sf = J2aeJ\f{i}^cn ^'^id A/'(i) is the set of links 
emanating from site i. The definitions of Fp, Pp are: 



(2) 



Pp EE /(5^, + 5^J • • /(5^3+5^J (3) 

where ai, q;2, aa, 014 are the links of plaquette p, enumer- 
ated clockwise. At a site with coordination number z, 
g{x) = {2x + z - 2)^/4. f{x) = 1 - x^. This model is 
closely related to the quantum dimer modeS (please see 
below) . 







FIG. 1. The action of the flip operator Fp on a typical 
plaquette. Notice that the total z-component of spin is gen- 
erally not conserved under such operation. The links with the 
up-spins are shown here as colored - this provides an alter- 
native graphical representation which will be exploited later 



The operator g (Sf) annihilates states (and only those 
states) which have S'f = —1, i.e. which have one and 
only one neighboring spin. The operator Fp 'flips' pla- 
quette p by flipping the four spins around it; an example 
of such flip is shown in Fig. The operator Pp is a pro- 
jection operator which annihilates all states except those 
in which up- and down-spins alternate around p - see 
Fig. H Plaquettes are assumed to have four sides, how- 
ever, they can be put together irregularly or can overlap, 
as the parallelogram-shaped plaquettes of the triangular 
lattice do. 

We will take Ji, J2, J3 > 0. 
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FIG. 2. The action of the projection operator Pp leaves 
plaquettes with the above shown spin configurations intact 
while annihilating any other type of plaquettes. A subse- 
quent application of the flip operator Fp to these plaquettes 
simply transforms them into each other, therefore they will 
be referred to as "flippable". 



III. THE LATTICE. 

Some care is required in the choice of lattice. As we 
will see below, the model (||) is tractable at J2 = J3. 
We would like to choose a lattice so that the ground 
state exhibits the key feature from which all of the inter- 
esting physics follows: fourfold ground state degeneracy 
on the torus even in the presence of local translational- 
symmetry-breaking fields. 

This can be accomplished if (1) the lattice does not al- 
low accidental symmetries which will increase the ground 
state degeneracy, a requirement which can usually be sat- 
isfied by taking a non-bipartite lattice and (2) the lattice 
has a unit cell which includes several plaquettes, so as 
to frustrate states in which the up-spins form an ordered 
crystal (i.e. a spin-density- wave) . There are many possi- 
ble lattices which satisfy these requirements. Our basic 
strategy for constructing these lattices is to take a Bravais 
lattice and introduce a periodic array of 'defects'. These 
defects pin a spin-density-wave state and make it non- 
degenerate, but they do not affect the fourfold degener- 
acy of the topologically-ordered state. We arrange these 
defects with a spacing which is incommensurate with the 
likely spin-density-wave states so that these states are 
frustrated and lifted in energy. Certain types of defects 
will also make it easier to satisfy (1). 

We will give two examples, T' and S'. T' is based 
on the triangular lattice (which, without defects, was ex- 
ploited in this context by Moessner and SondhiEZi, see be- 
low); the defects are missing sites, as depicted in Fig. ||. 
Even a single such missing site frustrates the staggered 
state - a special type of crystalline state with no "flip- 
pable" plaquettes - as depicted in Fig. ||. (A "flip- 
pable" plaquette is a plaquette which is not annihilated 
by Pp - see Fig. ||. It has alternating up- and down- 
spins whose direction can be reversed by application of 

/i.e.. Left to its own devices. 
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Fp will flip any plaquette, even the ones which are not 
'flippable'; Pp prevents this. On the triangular lattice, 
a plaquette is any primitive parallelogram.). In T' , an 
array of sites is missing, so that the lattice is given by 

|r R = ni(ax) + n2 ^^x + — ; "-1, "^2 ^ Omodfcj 

(4) 



where k is an arbitrary integer (this is just one such ex- 
ample, many other T'-type lattices can be constructed 
along these lines). Another possibility, S", is the square 
lattice iiL-which some of the plaquettes are split in two, as 
in Fig. EpI. This must be done so as to split sonic plaque- 
ttes horizontally and others vertically, in order to frus- 
trate staggered states aligned in both directions. These 
split plaquettes may be viewed as elementary dislocations 
in a perfect square lattice and thus they serve to model 
"real-life" defects. 
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FIG. 3. The mutilated triangular lattice, T' . The spins 
(only up-spins are shown) correspond to the maximally stag- 
gered configuration. In the presence of this type of lattice de- 
fects, there are strings of flippable plaquettes (shaded) which 
frustrate the true staggered state. 



IV. GROUND STATES AND TOPOLOGICAL 
ORDER. 

The first term in (^ requires each site to have one 
and only one neighboring up-spin, so that = —1, in 
which case it vanishes. On the square lattice, this leads 
to a magnetic moment which is half that of a fully po- 
larized state. In our model, this magnetic moment does 
not result from spontaneous ordering; the magnetization 
is actually fixed at a non-zero value. However, it is pos- 
sible to have a system which will spontaiieously undergo 
an Ising-like transition into such phasecS. We will call 
a state with both spin-charge separation and ferromag- 
netism F*, following the nomenclature of Ref. |^. The 
coexistence of conventional long-range order and quan- 
tum number fractionalization is familiar in ID and 2D: 
in polyacetylene, fractionalization coexists with CDW 
order; in easy-axis magnetic chains, with antiferromag- 
netism; at the v = 1/S quantum Hall plateau, charge 
e /3 quasiparticles can form a crystal and the topological 
order (and quantum Hall effect) will not be disrupted. 

The ground state of this model on T' or S' can be 
found exactly for J2 = J3. Every plaquette costs zero 
energy, so long as all flippable plajquettes are taken in 
the linear combination + Fp\^)c3. Since every spin 
configuration is obtainable frcipii every other one by the 
repeated application of FpPpESEJ, the ground state is the 
superposition with equal amplitudes of all possible con- 
figurations of spins satisfying 5f = —1. The ground state 
is annihilated by Hq. 
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Let us consider the crystalline states which compete 
with the topologiftally-ordered state. On the square and 
triangular latticed], the staggered state state does not 
mix with other states under FpPp since it is annihilated 
by this operator. It is a zero-energy ground state which 
is degenerate with the topologically ordered state. These 
two distinct ground states, have become degenerate at 
the first-order phase transition point J2 = Js- 

Fortunately, this is not the case on our lattices, as we 
now demonstrate. The staggered state has finite energy 
density at J2 = J3 on T' and S' since it is frustrated on 
these lattices. On T', there is no perfect staggered state. 
Consider a single defect (missing site). It is clearly im- 
possible to have a perfect staggered state in the presence 
of the defect. The closest that we can come to a per- 
fect staggered state (a which we will call a maximally 
staggered state) is either (a) a state which has a string 
of flippable plaquettes originating at the defect and ex- 
tending to infinity or (b) a state with one vertex which 
is frustrated by having = —2. If we take Ji — > 00, 
then only (a) is possible. When we introduce an array 
of defects, the strings of flippable plaquettes will orig- 
inate at one defect and terminate at another. Under 
the action of Fp, this state will mix with all of the oth- 
ers. Hence, a maximally staggered state will have energy 
density proportional to ^/p J2 in this limit (where p is the 
defect density). However, for J2 — > c», only (b) is pos- 
sible and the energy density of the maximally staggered 
state is proportional to p Ji . More generally, for Ji and 
J2 finite, the energy density of the maximally staggered 
state will be pJi/2 for small defect density p and will 
be proportional to y/pJ2 at large p. One might wonder 
whether there is some other crystalline state (e.g. one 
with a large unit cell) which has zero energy. However, 
if such a state contains 'flippable' plaquettes, it will mix 
under the action oi-^p with all of the other states with 
flippable plaquettea^y. Hence, such a state will cost finite 
energy. 

We can repeat the above analysis for S' . On S", we 
can frustrate one plaquette for each defect, with energy 
density pJ2- On the perfect square lattice, thaJ"* ground 
state is not fourfold degenerate. It is criticalEl and un- 
stable to a columnar state as J3 is decreasedl2^. This is 
not the case on T' or S", so we do not need to worry 
about the columnar state, either. 

The ground state may be visualized in the follow- 
ing way. Consider some reference configuration of spins 
which is annihilated by the first term of equation (Q) and 
color all of the links which have up-spins. Now take any 
other configuration which is also annihilated by the first 
term of equation (|l|) and do the same. By placing one 
graph on top of the other, and erasing all links at which 
both graphs coincide, we obtain a collection of loops on 
the lattice. If we visualize states in terms of their asso- 
ciated loop graphs, then the ground state is given by a 
superposition of different loop configurations. 

Since Hq conserves modulo 2 the winding numbers 
of these loops about either of the generators of the 



torus, there are 4 degenerate ground states on the torus, 
'0(ni,n2)j "-ij '^2 = 0, 1, with 0, 1 corresponding to even or 
odd winding numbers. By straightforward extension, the 
degeneracy on a genus g surface is 4^ . Although we have 
computed this degeneracy only at the special coupling 
constant value J2 = J3, we believe that it is robust over 
some range of parameters because it is characterized by 
an integer, 4. This integer cannot change as a result of 
infinitesimal perturbations, but only as a result of a per- 
turbation which is sufficiently strong that it moves the 
system across a phase transition at which this integer 
changes discontinuously. 

On the perfect square lattice, the directed winding 
number (not merely the winding number modulo 2) is 
conserved because the lattice is bipartite. As a result, 
there are Lx L sectors. S' is not bipactite, so it has only 
4 topologically-distinct ground states.E3 
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FIG. 4. A state of our model on the distorted square lattice 
with a spinon (centered about the site inside the dotted line), 
a holon (open circle), and a vison (a cross connected to a 
dashed line). The dotted-dashed line encloses one of the four 
defect plaquettes which has been split by additional sites. 



V. SPINONS AND HOLONS. 

The fourfold ground state degeneracy implies the ex- 
istence of fractionalized excitations, so long as it is un- 
related to translational and rotational symmetry, a con- 
dition which is satisfied as a result of our choice of lat- 
tice. To see this, imagine cutting open the torus along 
its second generator, thereby producing an annulus. The 
ground state "0(0,0) has non-zero projection on the ground 
state of the annulus because it will have some amplitude 
to have zero loops encircling the torus. '0(1, 0) does not, 
but it does have finite projection on a finite-energy ex- 
cited state of the annulus because it must have at least 
one loop circuiting the torus and it must have some am- 
plitude to have only one. This finite-energy excitation 
has, by construction, a spinon at the inner and outer 



4 



edges of the annulus, as in Laughlin's construction of 
charge e/3 quasiparticles in the fractional quantum Hall 
effects. The inner edge of the annulus can be shrunk 
and filled in since our discussion depends only on the 
topology, but not the geometry of the lattice. This con- 
struction can be done on a torus of any size, so the spinon 
at the boundary can be taken arbitrarily far away from 
the one in the interior with finite energy cost. 

This general argument can be substantiated in our 
model by a direct construction of spinon and holon exci- 
tations. Two spinous may be created by flipping a single 
up-spin into a down-spin. This changes by —1 and 
creates 2 sites with S'f = —2. These sites can be moved 
apart; each one carries = and costs energy Ji. 
A holon may be constructed by simply removing a spin 
from one of these spinon sites. This removes charge e 



from a neutral S*^ 



excitation. 



and spin 5^ = — ^ ^^^^^^ ^ ^^^^^^^^ ^ — ^ 
thereby producing a spinless charge e excitation. There 
will now be a site with S'f = —|, so the holon costs en- 
ergy Ji/4. In the loop picture, spinous and holons reside 
at the endpoints of broken loops. 



VI. VISONS. 



Consider now the operator 



(5) 



where Cp is any curve which starts at the center of plaque- 
tte p, connects it to the center of a neighboring plaquette 
p' , and continues in this manner through the centers of 
a sequence of neighboring plaquettes, running to infinity 
(or the boundary of the system). The product in equa- 
tion (|^) is over all links a which intersect Cp. Under the 
action of ^p, each loop configuration receives a — 1 if Cp 
has an odd number of intersections with colored links and 
1 if it has an even number of intersections with colored 
links. When a holon or spinon follows a trajectory encir- 
cling p, the intersection number must change by one, so 
$p creates a Z2 vortex, or 'visoncl, at plaquette p. 

The statistics of holons and spinous depend on the en- 
ergetics of the model: by binding to a vison, thejt-caii 
switch their statistics between bosonic and fermionicE3'E3. 
Our Hamiltonian does not allow holons or spinous to 
move: they are infinitely heavy. However, a small per- 
turbation will allow them to move and will give rise to 
the energetics which determines whether they bind with 
visons and, thereby, their statistics. 

The fourfold degeneracy of the ground state - or, in 
other words, the topological order - guarantees the ex- 
istence of a vison energy gap. To see this, consider the 
degenerate states ip(o,Q) i^(i,o) on the torus. Now imag- 
ine creating a vison pair at plaquette p, taking one vison 
around the second generator of the torus, and annihilat- 
ing the pair at p. This is equivalent to acting on our 
state with an operator similar to but with the curve 



Cp in equation replaced by a closed curve which passes 
through p and encircles the torus along its second gener- 
ator. This operator exchanges V'(o,o) i "0(1,0) • The ampli- 
tude for such a process is essentially the exponential of 
the Euclidean action required for such a virtual process 
to occur, ^ g-cLA„ jg length of the loop 

around the torus, A^, is the vison gap, and c is a con- 
stant. Hence, the energy splitting between states "0(0,0) 
and ^(1,0) is ~ e^"^^^". Since we know that this splitting 
vanishes in the thermodynamic limit, the vison gap Au 
must be finite. 

This conclusion is supported by a direct calculation. 
The creation of a vison at p takes the state \tp) + Fp\ip) 
into lip) — Fp\4>), with an energy cost A„. Since the vison 
creation operator <i>p commutes with all of the terms in 
equation (!]) except for the J2 term at plaquette p, with 
which it anticommutes, a state with one vison, has 
excitation energy: 



H \%) = 2 J2 



Pp\0) 



(6) 



Hence, the vison gap is equal to 2J2 multiplied by the 
density of flippable plaquettes. This may be computed at 
J2 = J3 by the Grassmann techniques discussed below. 
In an integrable model which we discuss below, exact 
vison eigenstates and energy eigenvalues may be found. 

From equation (^), we see that the vison gap will be 
non-vanishing whenever (Oj-PpPplO) 7^ 0, i.e. whenever 
the spins fluctuate in the ground state, as they gener- 
ically do in our model, even outside the topologically 
ordered phase. However, this is not particularly conse- 
quential. Consider the analogous situation in a super- 
fluid: it is possible to define a vortex energy above the 
transition (e.g. the Kosterlitz-Thouless transition) which 
varies smoothly across the transition. However, this en- 
ergy is only meaningful in the superfluid state (or, per- 
haps, near it). Similarly, the vison gap becomes mean- 
ingful in the topologically ordered phase. Outside this 
phase, $p is merely an operator which creates some com- 
plicated gapped excitation. 

A vison gap is necessary for topological order; a gap 
in the rest of the spectrum is not. However, the equal- 
time spin-spin correlation function in the ground state 
is exponentially decaying, as may be seen from an exact 
mapping between the ground state_Df our model and the 
field theory of free lattice fermionscJ, according to which 
it is a square root of the eight-fermion correlator. This 
decays exponentially with distance since these fermions, 
though massless on a regular square lattice, acquire mass 
in the presence of lattice distortions such as that shown 
in Fig. I (the details of this calculation will be published 
elsewhere). Hence, it is natural to conclude that there is 
an energy gap to all excited states; this may be argued 
via the single-mode approximation. This computation 
has been carried-put on the triangular lattice by Moess- 
ner and Sondh£j, who found that the single-mode ap- 
proximation suggests that the system is indeed gapped. 
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VII. A FAMILY OF MODELS. 



IX. QUANTUM DIMER MODELS. 



We can gain further insight into the spectrum of our 
model by generalizing it to the following family of models: 

i 

-J'^\Fv (l-Pp)-(l-Pp)] (7) 
p 

where y(x) = |a; {x^ — 4) (a: + 1) (for z = 4) and J2 — J3 
in Hq. As J{, J2 are iiicreased, two things occur: loop 
crossings are allowedEj and plaquette flips of 'unflip- 
pable' plaquettes are allowed. The topological order is 
preserved because the winding number modulo 2 is con- 
served. In the extreme limit J[ = Ji, J2 = J2, there 
are eight equally likely configurations at each site (cor- 
responding to those of the eight-vertex model) and the 
model is now integrable because [g (S'f ) -I- y (Sf)] and Fp 
commute among themselves and, hence, with the Hamil- 
tonian, for all i, p. Hence, we can simultaneously diago- 
nalize all of these operators. They have eigenvalues and 
1, respectively, in the ground state; a g (S'f) -I- y {Sf) = 1 
eigenvalue is a quasiparticle excitation at i and Fp = =J 
is a vison at p. This model is equivalent to Kitaev'sciJ. 
The ground state of this integrable model has the same 
topological order (fourfold degeneracy) as that of (|l|), but 
its quasiparticle excitations do not carry spin since it is 
not conserved. Since there is no projection operator Pp 
in H, the vison energy is exactly 2J2- Crystalline states 
have energy J2/2 per plaquette above the ground state. 

VIII. FIELD THEORETIC DESCRIPTION. 

The configurations allowed in the ground state of the 
integrable model are described by closed loops - in other 
words, by the configurations of the Ising model on the 
dual lattice (equivalent to the eight-vertex model). The 
dynamics of the plaquette flip operator is the same as 
that of a transverse field in the Ising model. Hence, 
the integrable model is equivalent at low energies to the 
transverse field Ising model which, in turn, is dual to a 
Z2 gauge theory. Since the topological order associated 
with Hq is the same as that of the integrable model, it, 
too, is described by a Z2 gauge theoryllj, as proposed in 
refs |0||. 

Note that Sf is conserved for all i in our model (|^). 
Hence, there is an independent U(l) symmetry at each 
site of the lattice. However, only time-independent trans- 
formations leave the Lagrangian invariant. Hence, this is 
an ordinary (but large) symmetry group; it is not a gauge 
symmetry, which must allow time-dependent transforma- 
tions. All of the degrees of freedom in (|]) are physical. 
This is similar to the symmetry of a set of noninteracting 
spins in a magnetic field, H = — X);-^ ' which also has 
an independent U(l) at each site of the lattice. 



Ouij^uindiiL(^ can be mapped to the quantum dimer 
modeH'E3E3o, in which it is assumed that that there are 
spins located at the sites of a lattice and that each spin 
forms a singlet dimer with one of its nearest neighbors. 
In (^, an up-spin link corresponds to a dimer; a down- 
spin link to the absence of a dimer; spinous, to empty 
sites (which are holons in the dimer model). Then the 
first term in equation (^ requires each spin to form a 
dimer with exactly one of its neighbors. The J2 and J3 
terms are precisely the dimer kinetic and potential ener- 
gies of Ref. |2^. Our F* StMe of (|l|) is simply the resonat- 
ing valence bond (RVB)ql3 ground state of the quantum 
dimer modeLon the same lattice. Recently, Moessner 
and Sondh£j gave compelling evidence that the trian- 
gular lattice quantum dimer model has an RVB ground 
state over a substantial range of parameters terminating 
at a first-order phase transition at J2 = J3 into the stag- 
gered state. According to our arguments, the RVB state 
is the unique, exact ground state at J2 — J3 on T' . 

X. SU(2) SYMMETRIC MODELS. 

There is no reason to believe that an SU(2) symmet- 
ric magnet cannot be topologically-ordered. If we wish 
to apply the preceding results, then the quantum dimer 
model can be the starting point for the discussicpnpf||aii 
SU(2) symmetric topologically-ordered magnetcS'EZlEa. 
There is some numerical evidenced that a model of 
spins on the triangular lattice with strong four-spin ex- 
changes - reminiscent of the plaquette flip operator - has 
a topologically-ordered ground state with fourfold degen- 
eracy. 

XI. SUMMARY. 

We have demonstrated that quantum number frac- 
tionalization and topological order are an eminently rea- 
sonable possibility for two-dimensional magnets by con- 
structing a microscopic model which exhibits these phe- 
nomena. Our result links two problems of great interest: 
quantum number fractionalization in 2D quantum mag- 
nets and the investigation of physical systems which are 
suitable platforms for fault-tolerant quantum computa- 
tion (this link was also noted in Ref. A possible 
nexus with ideas about superconductivity - either in the 
cuprates or elsewhere - leads to potentially fruitful av- 
enues for further research in both areas. 
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